Metamagnetism of antiferromagnetic XXZ quantum spin chains 



Toru Sakai 

Faculty of Science, Hvmejt Institute of Technology, Ako, Hyogo 678-1297, Japan 

Minoru Takahashi 

Institute for Solid State Physics, University of Tokyo, Roppongi, Minato-ku, Tokyo 106-8666, Japan 

(March 1999) 

The magnetization process of the one-dimensional antiferromagnetic Heisenberg model with the 
Ising-like anisotropic exchange interaction is studied by the exact diagonalization technique. It 
results in the evidence of the first-order spin flop transition with a finite magnetization jump in the 
Neel ordered phase for S > 1. It implies that the 3 = \ chain is an exceptional case where the 
metamagnetic transition becomes second-order due to large quantum fluctuations. 
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Metamagnetism is one of interesting topics in the field 
of the magnetism. The spin flopping of the antiferromag- 
nets with the Ising-like anisotropic exchange interaction 
is a simple mechanism of the field-induced metamagnetic 
transition. § In the process the spins abruptly changes 
directions from parallel to perpendicular with respect to 
the easy axis of the sublattice magnetization at some 
critical magnetic field. In the classical spin systems it is 
easily understood that the transition is first-order and ac- 
companied by a finite magnetization jump independently 
of the lattice dimension. In the quantum systems, how- 
ever, the process is not so trivial because the sublattice 
magnetization shrinks or vanishes due to quantum fluc- 
tuation even at zero temperature. Such quantum effect is 
expected to be larger in lower dimensions. In one dimen- 
sion (ID) the 5 = | Ising-like XXZ model was proved 
by the Bethe ansatz exact solution to exhibit a second- 
order metamagnetic transition at some critical magnetic 
field. In contrast the recent numerical analysis [|3| sug- 
gested that the correponding transition is first-order for 
the same model in two and three dimensions (2D, 3D), 
while the magnetization jump has a quantum spin reduc- 
tion. It implies that the quantum fluctuation plays such 
an important role as to change the order of the phase 
transition in low dimension. As far as we restrict us on 
the 5 = I model, the critical dimension seems to lie be- 
tween one and two. On the other hand, the strength of 
the quantization also depends on the spin value 5. Since 
the quantum effect is smaller for larger 5, the second- 
order transition descovered for the 5 = \ chain is not 
necessarily a common feature for all the values of 5. At 
least the transition is first-order in the infinite 5 limit. 
Thus it is important to investigate the critical value 5 C 
between the second- and first-order transitions, for un- 
derstanding the quantum effect on the metamagnetism 
in ID. In this paper, we study the magnetization process 
of the 5 = 1 and | XXZ spin chains, particularly on the 
possibility of the first-order spin flopping transition. 

The magnetization process of the ID XXZ model is 
described by the Hamiltonian 



TL = 7io + TLz, 

3 

Hz= —H Sj, 

j 

under the periodic boundary condition. We restrict us 
on the Ising-like anisotropy < 7 < 1 to consider the 
metamagnetic transition from the Neel ordered ground 
state, induced by the external magnetic field along the 
easy axis of the sublattice magnetization. In the 5 — i 
chain a second-order transition occurs at some critical 
field H c for < 7 < 1. The asymptotic behavior has the 
form 

m~{H-H c ) 1/2 . (2) 

The magnetization curves in the thermodynamic limit 
for 7 =0.3 and 0.5 calculated by the Bethe ansatz are 
shown as solid lines in Fig. |], where H s is the satu- 
ration field (H s = 25(1 + 7)). On the other hand, in 
the limit of 5 — > 00 we know that the spin flop oc- 
curs at H c = H s [(l - 7)/(l + 7)] 1 / 2 . This is the first- 
order transition and the magnetization jumps from to 
5[(1 — 7)/(l + 7)] 1 / 2 , as shown in Fig. g. In order to in- 
vestigate the corresponding transition for finite 5 larger 
than 1/2, we perform the exact diagonalization technique 
applied to the finite clusters with the system size L, be- 
cause there is no exact solution available here. Using 
the Lanczos algorithm, we calculate the lowest eigen- 
value of Hq in the subspace with ^ ■ 5| = M for the 
L-site system, defined as E(L, M). To derive the magne- 
tization curve, the magnetic field H to bring about the 
magnetization m = ^ and ml = m + | is estimated by 
[E(L, M+1)-E(L, M-l)]/2 and E(L, M+1)-E(L, M), 
respectively. To demonstrate the validity of the method, 
we show the results of the 5 = \ chain by the diagonal- 
ization for L = 20 and 24 as circles, together with the 
Bethe ansatz solution in Fig. |l|. The bulk property is re- 
vealed to be well understood even in such small systems. 

In the small magnetization limit in the Neel ordered 
phase, only the state with M — 2Sn (n = 0, 1,2, ■■•) 



1 



can contribute the ground-state magnetization process 
of the spin-.? chain. Because states where every site has 
S? = +S or — S are stabilized better than other states. 
To explain the reason, we consider the S = 1 chain, for 
example. In the excited state produced by changing Sj 
from —1 to in the Neel state (• • • IT tl 4 • ■)- the mo- 
tion of the object to the adjacent site (• • • |tt | • • •) 
will lead to some energy loss in the diagonal element of 
the exchange Hamiltonian, while it will gain some en- 
ergy in the off-diagonal element. In contrast the object f 
doped in the Neel order (• • • ITtTIT ' ") can yield the en- 
ergy gain due to the motion along the chain with no loss 
in the diagonal element (••• ITTITT •'•)■ Thus smaller 
magnetic field can excite J, to f , rather than to 0. It im- 
plies that states with +1 or —1 at every site are more 
favorable than states including 0, and the magnetization 
process starting from the Neel phase consists of the states 
with M = 2n in the small m limit. The argument is also 
generalized for arbitrary S and the state with M = 2Sn 
is stabilized. Thus in the small magnetization region an 
oscillation with the period of 2S in the M dependence 
of the energy of the spin-5 chain in the Ising phase. 
To avoid the oscillation, we will take only states with 
M = 2Sn into account in the following analysis, except 
for the determination of the saturation field H s . 




FIG. 2. Magnetization curves in the classical limit. 
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FIG. 1. Magnetization curve for S = | obtained by the 
Bethe ansatz (solid lines), and exact diagonalization of the 
finite chains with L = 20 and 24 (symbols). 



FIG. 3. Schematic curve of e(m) in the case when the 
first-order metamagnetic transition occurs. m a is the magne- 
tization jump and it is determined by eo(m) — 0. 

In general, a first-order metamagnetic transition oc- 
curs, when the ground state energy per site e(m) for the 
Hamiltonian Ho (We choose the origin to set e(0) = 
here.) has such m-dependence as shown in Fig. [3|, where 
the region with e"(m) < lies from the origin to some 
magnetization. We define eo (m) as the intersection of the 
vertical axis and the tangent line of the curve e(m), as 
shown in Fig. [| The magnetization jump m s is deter- 
mined as the solution of eo(m) = 0. It implies that the 
sign of eo(m) changes at m s and the magnetization with 
eo(m) > does not appear in the ground-state magneti- 
zation process. In order to examine the possibility of the 
magnetization jump in the spin-S* chains, we investigate 
the behavior of eo(m) based on the finite cluster calcula- 



2 



tion. We assume the size dependence of the ground state 
energy per site obeys the relation 



e(L,M) = y[E(L,M)-E(L,0)] 



e (m) + 0(-l), 



(3) 



for L — > oo with fixed m — Although the relation is 
predicted by the conformal held theory Q for massless 
state, the following argument will be valid even in mas- 
sive cases, as far as the size correction does not exceed 
O(-Jj) We also dehnc h(L,M) as 



h(L, M) 



1 

25 

e (m) 



[E(L, M + 25) 



E(L,M -2S)} 



(4) 



where the size dependence is derived from (||) . The quan- 
tity h(L, M) should converge to H for L — > oo in the 
normal magnetization process. Since eo(m) is given by 
eo(m) = e(m) — e'(m)m, it can be estimated by the form 



L 



[e(L, M) - h(L, M)M] ~ e (m) + 0( 



L 2 ' 



(5) 



Since various system sizes with fixed m — are avail- 
able only for few values of m because of the restriction 
A/ = 2Sn, we neglect the size correction in (||) in the fol- 
lowing analysis. For 5 = i the finite cluster calculation 
up to L = 26 indicates eo(m) < for < m < 1 in all 
the region < 7 < 1. It means no first-order transition, 
which is consistent wiht the result from the Bethe ansatz. 
For S = 1, however, the calculation up to L = 20 detects 
the region with eo(m) > 0, which suggests the first-order 
metamagnetic transition. When we vary 7 with fixed M , 
a point with eo{m) = can be found for small m. The 
points for various values of M available for L = 18 and 
20 are plotted on the m/S-j plane as circles in Fig. ^. 
Since ea(m) is positive under the points, while negative 
over them in the plane, the points stand for the magne- 
tization jump m s for corresponding 7. In the analysis 
of 20-site cluster, eo(m) is always positive for M > 12. 
Thus we think that m = 0.6 is an upper bound of m s in 
the limit 7^0 denoted as to„o, which is shown as x in 
Fig. ||. m s0 does not correspond to the value of the Ising 
model m s /S = 1, which implies that there exists a spin 
reduction due to the quantum fluctuation, as was found 
in the 2D and 3D S = 5 model. || On the other hand, 
the 5 = 1 XX Z chain has the Haldane phase P3] , where 
the second-order transition described by (0) was revealed 
to occur in the magnetization process |P- f7|, in the large 
7 region. Consulting some previous works iflq [Tl|| , we put 
the boundary between the Haldane and Nccl phases at 
7 C = 0.84, as shown in Fig. ^. In the present work, by 
the analysis of several excitation gaps corresponding to 
the soft modes of the Luttinger liquid, we checked that 
the magnetic state for m s < m < 1 consists of a single 
gapless phase characterized by the dominant spin corre- 
lation function (SqS~) ~ (— l) r r~' n , in contrast to the 



magnetization process in the presence of the Ising-like 
single-ion anisotropy described by DJ2j(Sj) 2 [D < 0), 
where there appears another massless phase ||?] with 
(5q5^) ~ cos(2kFr)r~' n . Thus the line of m s is expected 
to continue to the unique phase boundary j c at m = 0. 
The solid line for 5 = 1 in Fig. || is the result from the 
polynomial fitting to the points with eo(ro) = (L = 18 
and 20) and the boundary 7 C , also adding m s o = 0.6 be- 
cause the fitted line excluding it would exceed the upper 
bound. The points with eo(m) = for 5 = | are also 
plotted as squares in Fig. [|, where we use the results of 
L = 12 and 14. Assuming the boundary of the Neel phase 
lies at 7 = 1, the curve fitted to the points is shown as a 
solid line in Fig. ^. In this time m s o is estimated from 
the fitted curve as m s o = 0.70±0.04, because it is smaller 
than the upper bound determined by m for which eo(m) 
is always negative. The spin reduction in the limit 7^0 
is smaller than that for 5=1. It is consistent with the 
reasonable assumption that the quantum fluctuation is 
smaller for larger 5. The line of the magnetization jump 
in the classical limit m s /S — [(1 — 7)/(l + 7)] 1 ^ 2 is also 
plotted as a dashed curve in Fig. ^. These curves sug- 
gest that the magnetization jump increases with the spin 
value increasing from 5 = 1 towards the classical limit. 
Therefore we conclude that the ID Ising-like XX Z model 
(0 < 7 < 1) exhibits the first-order spin flopping transi- 
tion for 5 > 1, except for the Haldane phase. It implies 
that the 5 = | chain is an exceptional case where the 
metamagnetic transition from the Neel phase is second 
order, that is, the critical spin value 5, 
and 1. 
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FIG. 4. Magnetization jump m s /S for S = 1 and | (solid 
lines) obtained by the polynomial fitting to the points with 
6(0) = for the finite systems (symbols). For S = 1 7 C is the 
boundary between the Neel and Haldane phases, and x is the 
upper bound of m S Q. m s /S in the classical limit is a dashed 
line. 
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FIG. 5. Magnetization curves for S = 1. 
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FIG. 6. Magnetization curves for S 1 



In order to convince of the conclusion, we show the 
ground-state magnetization curves for 5 = 1 and | 
based on the present finite cluster calculation. We plot 
h(L,M) as the field H for the magnetization m = 
and h'(L,M) = [E(L,M + 25) - E(L, M)]/(25) for 
m = (M + S)/L, using the value renormalized by the 
saturation field H s . The results of the largest two sys- 
tem sizes for 5 = 1 and J (7 = 0.3 and 0.5) are shown as 
symbols in Figs. || and [| respectively. They all indicate 
the evidence of the first-order transition, in contrast to 
the 5 = I in Fig. [I]. We determine the magnetization 



jump m s from the fitted lines in Fig. 4, and put the fit- 
ted curves for m s < m in Figs. |^ and 6, to complete the 
magnetization curves. The behaviors of the symbols in 
Figs. H and || suggest that the size correction is so small 
that the magnetization curves given by the solid lines are 
expected to be of the bulk systems. 

The present analysis indicated a numerical evidence of 
the first-order metamagnetic transition in the ID XX Z 
model for 5 > 1, based on the finite-cluster calculation. 
In general the finite system has larger quantum fluctu- 
ation, which tends to suppress the magnetization jump, 
than the infinite system. Thus the present conclusion 
on the finite jump is expected to be valid even in the 
thermodynamic limit. 

In summary, the exact diagonalization study on the fi- 
nite cluster suggested that the ID Ising-like antiferromag- 
netic XX Z model exhibits the first-order spin flopping 
transition with a finite magnetization jump for 5 > 1, 
except for the Haldane phase. It implies that only for 
5=5 the quantum fluctuation is extremely large and 
makes the corresponding transition second-order. 

The numerical computation was done using the facil- 
ity of the Supercomputer Center, Institute for Solid State 
Physics, University of Tokyo. 
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